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The sensitivity of an experiment to detect the Majorana neutrino mass via neutrinoless double beta decay 
(Ov 6B) strongly depends on the rate of background events that can mimic this decay. One major source of this 
background is the radioactive emissions from the laboratory environment. In our study, we focused on assessing 
the background contributions from environmental gamma rays, neutrons, and underground muons to the Jinping 
bolometric demonstration experiment. This experiment uses an array of lithium molybdate crystal bolometers 
to probe the potential 0186 decay of the °°Mo isotope at the China Jinping Underground Laboratory. We also 
evaluated the shielding effectiveness of the experimental setup through an attenuation study. Our simulations 
indicate that the combined background from environmental gamma rays, neutrons, and muons in the relevant 
100Mo 0v 86 Q-value region can be reduced to approximately 0.003 cts/kg/keV/yr. 
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I. INTRODUCTION 


Unveiling the properties of neutrinos is fundamental, and 
searching for Neutrinoless double beta (0v 86) decay stands 
as a pivotal approach to this goal. Observing the 0v 6 decay 
experimentally can directly substantiate the Majorana nature 
of neutrinos, implying that neutrinos are their own antipar- 
ticles. This would open new frontiers in physics, extending 
beyond the standard framework of particle physics [1]. 

Numerous experimental techniques have been employed 
to detect the rare OvG( events. While direct evidence of 
Ov 6 remains elusive, strict limits on the half-life of var- 
ious double beta decay nuclides have been established [2— 
5]. Among the tools used, crystal bolometers, which operate 
at ultra-low temperatures, stand out. They are preferred for 
Ov searches due to their high energy resolution, efficient 
detection, and stable operation [6]. These cryogenic bolome- 
ters employ crystals embedded with double-beta decay iso- 
topes to detect the 0v8 decay. The forthcoming CUPID 
experiment (CUORE with Particle Identification), integrating 
a high-performance heat-light dual readout to discern alpha 
particle backgrounds, is projected to enhance the sensitivity 
of Ov G6 detection in the future [7-9]. 

A critical challenge in 018 detection is the potential in- 
terference from external radioactivity, as it can significantly 
hamper the experiment’s potential. Efforts have thus been 
made to understand and mitigate sources of background in- 
terference. Previous research has identified multiple culprits 
for misleading counts in the 0v{6 search energy region of 
interest (ROI), including cosmogenic activation, natural ra- 
dioactive contamination from materials, and environmental 
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radioactive emissions in the laboratory [10-14]. Given that 
cosmogenic activation at sea level notably influences the de- 
tectors, experiments are usually conducted in deep mines 
or under mountains; these locations, shielded by the earth’s 
crust or mountain rock, offer protection from cosmic radia- 
tion [15, 16]. Besides residual cosmogenic background, envi- 
ronmental radioactive interference remains crucial for under- 
ground setups. Main sources include gamma rays from rocks, 
neutrons from surrounding materials, and events induced by 
underground muons [17, 18]. Passive shields are typically 
deployed to counter gamma and neutron backgrounds, while 
muon vetoes are used to effectively tag and exclude under- 
ground muon events and their induced activities. 

Given that background levels are pivotal for experimen- 
tal sensitivity, thorough studies on environmental background 
shielding and subsequent evaluations of any remaining back- 
ground are crucial during the design phase of the detector 
system. In this paper, leveraging a Geant4-based Monte 
Carlo simulation, we delve into the background contribu- 
tions—due to environmental gamma rays, neutrons, and cos- 
mic muons—to the Jinping CUPID demonstration experi- 
ment. This experiment, featuring a 10 kg cryogenic crys- 
tal bolometer at the China Jinping Underground Laboratory 
(CJPL), aims to showcase the advanced technologies intended 
for the next-generation bolometric experiment targeting the 
OvGB decay of 1°°Mo. This paper is structured as follows: 
the simulation setup and methods are discussed in Sect. IL. 
Findings related to particle attenuation and background are 
elaborated in Sect. III. Conclusions and discussions can be 
found in Sect. IV 


Il. SIMULATION SETUP 


We utilized the Geant4 simulation toolkit (version 4.10.5) 
to model radiation shielding and assess background contri- 
butions [19]. Geant4 is renowned for offering an extensive 
array of particle-material interaction cross sections, allowing 
precise characterization of particle production and transport 


properties [20-22]. For our study, we employed the "Shield- 
ing" physics list, which is tailored for simulations associ- 
ated with deep underground rare event experiments. This 
physics list incorporates a series of interaction models, each 
detailing specific physical processes. For neutrons with ener- 
gies up to 20 MeV, we leveraged high-precision models and 
cross-sections to depict processes such as elastic and inelas- 
tic scattering, capture, and fission. To study muon-induced 
backgrounds, we also applied the "QGSP_BERT" physics 
list as a supplementary check, given its suitability for cos- 
mic ray-related processes. Event generation for gamma rays 
and neutrons was based on the in-situ spectral measurements 
taken in the CJPL experimental hall, as depicted in Fig. 1. 
We assumed isotropic radiation for both gamma and neutron 
sources. The initiation points for these environmental back- 
ground sources were set in line with the hall’s dimensions (14 
m in height and 12 m in width). Moreover, the muon sim- 
ulation factored in experimentally measured zenith and az- 
imuthal angle distributions. 
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Fig. 1. (Color online) Energy spectra of CJPL environmental gamma 
(left) and neutron (right). The sampled Monte Carlo spectra (solid 
lines) are shown in comparison with the results obtained by the in- 
situ measurements (dashed lines) [23, 28]. 
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Fig. 2. Schematic view of the cryogenic system, shielding system 
(a), and crystal detector array (b) implemented in the Monte Carlo 
simulation. 


Fig. 2 shows a schematic view of the experimental setup, 
including the detector crystal array, background shields sur- 
rounding the detector, and cooling system used to main- 
tain low temperatures. The configurations of the detector 
and shielding systems are implemented in the simulation 


Table 1. Key elements of the experimental setup. 


Element Material Mass (kg) Surface 
(cm?) 
Crystal array LMO 10.0 4374.0 
Thermal shield Cu-OFE 61.7 23850.9 
Top Cu shield Cu-OFE 209.7 4736.9 
Lateral Cu shield Cu-OF 2990.4 60880.0 
Lateral PE shield PE 1343.9 114300.0 


code. The CUPID-CJPL demonstrator comprised an array 
of thirty-six lithium molybdate (Li2MoO4) or LMO cryo- 
genic bolometers, having a cumulative crystal mass of ap- 
proximately 10 kg. The LMO crystals not only serve as the 
source for the 0v 88 decay but also function as heat absorbers. 
This array is organized into nine layers, with each layer con- 
taining four bolometers. Each of these bolometers is cooled 
to an ultra-low temperature of about 10 mK using a dilution 
refrigerator. The individual bolometer module integrates a 
4.5x4.5x4.5 cm? LMO cubic crystal and is equipped with 
a heat-light dual-phase readout. For these modules, 280 g 
LMO crystals were specifically chosen. The cooling appara- 
tus, or the cryostat, consists of three concentric vessels dedi- 
cated to thermal shielding. To complement these thermal bar- 
riers, radiation shields were strategically positioned both in- 
side and surrounding the cryostat to counteract radiation em- 
anating from the cryogenic system. Within this arrangement, 
a substantial copper shield (having a thickness of 12 cm and 
a diameter of 50 cm) was situated above the detector tower to 
mitigate external radioactive interference. To further fortify 
the system, passive shielding materials were employed. This 
includes a 12 cm thick layer of copper, designed to diminish 
environmental gamma rays, and a 20 cm thick layer of 5%- 
borated polyethylene (PE) tailored to intercept environmental 
neutrons. The salient features of this experimental setup are 
detailed in Table 1. 


MI. RESULTS AND DISCUSSION 
A. Attenuation study 


Numerous experimental studies have established that ma- 
terials with high atomic mass numbers, such as lead and cop- 
per, excel at attenuating gamma rays. Notably, copper typi- 
cally possesses lower inherent radioactivity than lead. As a 
result, high-purity oxygen-free copper has emerged as a pre- 
ferred material for crafting shielding systems in sensitive un- 
derground experiments related to rare events. 

For the Jinping bolometric demonstration, the design of the 
background shielding system was informed by a Monte Carlo 
simulation. Utilizing the Geant4 toolkit, we simulated parti- 
cle interactions with materials to quantitatively assess the at- 
tenuation capacities of the shielding against gamma rays and 
neutrons. 

Fig. 3 shows the gamma-ray absorption fraction as a func- 
tion of the gamma energy for different copper shielding thick- 
nesses. The incident gamma energy was sampled based on 
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Fig. 3. (Color online) Simulated results of the gamma absorption for 
varying thicknesses of the copper brick. 
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Fig. 4. (Color online) Simulated results of the attenuation coefficient 
as a function of copper shielding thickness for gamma and borated 
polyethylene thickness for neutron. The input energy spectra for 
gamma and neutron are shown in Fig. 1. 


the experimentally measured CJPL in situ gamma spectra. 
1. Our results indicate that the gamma absorption at differ- 
ent energies improves significantly by increasing the thick- 
ness of the material. In addition to high-energy environmen- 
tal gamma rays, fast neutrons with energies above 1 MeV are 
serious background sources for rare underground event ex- 
periments. Borated polyethylene (BPE) is an effective ma- 
terial for fast neutron shielding [24, 25]. The results of the 
neutron-shielding performance are shown in Fig. 4 compared 
with the gamma results. The black and blue curves represent 
the gamma and neutron absorptions of copper and 5% BPE, 
respectively. It was observed that the neutron shielding was 
improved by increasing the thickness of BPE. A reduction in 
gamma and neutron fluxes by a factor of ~10 can be real- 
ized by deploying a 12-cm thick copper shield and an 8-cm 
5% BPE shield, respectively. Thus, for the CJPL bolometric 
experiment, with 12-cm thick copper and 20-cm thick 5%- 
borated polyethylene external shields, we expect a reduction 
exceeding 10-fold in the environmental gamma and neutron 


fluxes within the detector system. 


B. Environmental gamma induced background 


Different background sources have the potential to simu- 
late a Ov GP signal within the decay-energy ROI of 0v 8. 
A predominant source of such radiation is environmental 
gamma rays. For underground experiments, these gamma 
rays largely stem from the decay of radioactive nuclei and 
processes induced by high-energy neutrons. Moreover, the 
lab’s infrastructure and its experimental apparatus can in- 
troduce additional gamma background. In our research, 
we delved into evaluating the contribution of environmental 
gamma background to a !°°Mo 0v 86 detection using a bolo- 
metric method. 

The gamma radioactivity of the CJPL environment has 
been extensively profiled in prior studies through in-situ 
gamma-ray measurements employing a high-purity germa- 
nium spectrometer [26-30]. The naturally occurring gamma 
spectrum spans up to 3 MeV; events exceeding this energy are 
typically the result of interactions between neutrons, muons, 
and rocks. Given that the isotope !°°Mo boasts a relatively 
high Qgg value (approximately 3034.40(17)keV), we des- 
ignated the ROI to range from 3000-3060 keV. This deci- 
sion factored in the energy resolution of the crystal bolometer 
when assessing the background. The approach to parameter- 
izing the energy resolution is as follows: 


o(E) = 0.49 + (0.058 x VE)?, (1) 


where E denotes the measured energy in keV based on the 
CUPID-Mo study [31]. 

The gamma spectrum measured using a high-purity germa- 
nium (HPGe) detector in CJPL Phase II was used as the input 
for our Geant4 simulation [26]. In this study, given that we 
focused on the °°Mo 08 ROI, only high-energy gamma- 
rays (EL, >3000 keV) were considered. The integral flux of 
the high-energy gamma rays was calculated as 2.51 + 0.07 
x 1076 em~?s~! based on the experimentally measured re- 
sults. 

In bolometric experiments, the array structure of the detec- 
tor system allows us to apply the anti-coincidence cut crite- 
rion. This helps in rejecting high-energy external background 
events that pass through several crystal modules. Fig. 5 il- 
lustrates the concept of “multiplicity”. An event that de- 
posits energy in a single crystal module is defined as hav- 
ing “multiplicity = 1”. Conversely, if the energy is de- 
posited across multiple crystal modules, it’s classified as a 
“multiplicity > 1” event. By selecting the “multiplicity = 
1” events, we can significantly filter out high-energy exter- 
nal gamma rays, as well as neutrons and muons, that traverse 
multiple crystals. 

Fig. 6 presents the simulated results for the gamma- 
induced background. The results for different cutting criteria 
are shown for comparison. It is observed that by requiring 
“multiplicity = 1”, the environmental gammas contribute 


Multiplicity = 1 


LMO 


Multiplicity = 2 


—o-f 


Fig. 5. (Color online) Definition of the multiplicity. 
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Fig. 6. (Color online)Simulated results of the gamma-induced back- 
ground. 


approximately 1.3 x 1073 cts/keV/kg/yr background in the 
100Mo 0v 68 ROI (3000-3060 keV). The uncertainty of the 
result is estimated by considering the statistical uncertain- 
ties of the flux measurements and selection of PhysicsList 
(“Shielding” vs “QGSP_BERT_HP”). The relative difference 
was found to be within 10%. 


C. Environmental neutron induced background 


High-energy neutrons present another significant back- 
ground for Ov 8p experiments. In deep underground labo- 
ratories, neutron fluxes below 10 MeV primarily arise from 
spontaneous fission (such as from U-238) and (a,n) reac- 
tions, which involve interactions between alphas from natu- 
ral radioactive sources and light target nuclei within the rock. 
Beyond production by the inherent radioactivity of materi- 
als, neutrons can also be produced by cosmic-ray muons and 
their secondary particles through nuclear reactions [32, 33]. 
The energy of these muon-induced neutrons can reach up to 
several gigaelectronvolts. 

Systematic measurements of the fluxes of low-energy ther- 
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Fig. 7. (Color online) Simulated results of the neutron-induced back- 
ground. 


mal neutrons and high-energy fast neutrons have been per- 
formed in deep underground laboratories [23, 34—36]. In 
our simulation study, we evaluated the neutron-induced back- 
ground based on the in situ fast neutron spectrum measured 
at the CJPL Hall-A [37]. Monte Carlo neutrons with kinetic 
energies exceeding 1 MeV were sampled based on the ex- 
perimental flux spectrum results. The total flux was approxi- 
mately (1.51 + 0.03) x 1077 em~?s7! for E,>1 MeV with- 
out the polyethylene shielding. The results for the residual 
neutron-induced background with the copper and PE shield- 
ing setups are shown in Fig. 7. We observed that the en- 
vironmental neutrons contribute approximately 1.6 x 107° 
cts/keV/kg/yr in the °°Mo 0v 66 ROI with “multiplicity = 
1” cut applied. It should be noted that there are currently no 
published experimental data on CJPL environmental neutron 
flux exceeding 20 MeV. Thus, high-energy neutrons (E, >20 
MeV) were not considered in our simulations. Neverthe- 
less, according to existing simulation studies performed for 
LNGS, the integral flux exceeding 10 MeV (mainly from 
muon-induced processes) is three orders of magnitude lower 
than that below 10 MeV (mainly from radioactivity-induced 
processes) [38, 41]. For the CJPL, owing to the larger depth, 
the flux of cosmic-muon-induced neutrons was even lower. 
The relative background contribution from EF, >20 MeV high 
energy neutron was expected to be low. 


D. Underground cosmic muon induced background 


Cosmogenic backgrounds play a pivotal role in rare-event 
search experiments [39]. At Earth’s surface, cosmic rays pre- 
dominantly consist of high-energy muons [40]. The signifi- 
cant rock overburden of CJPL (equivalent to 6700 m of water) 
dramatically reduces the cosmic muon flux underground—by 
about eight orders of magnitude compared to the flux at sea 
level. As a result, the cosmogenic background doesn’t pose a 
significant threat to the Jinping bolometric experiments. Nev- 
ertheless, for the heightened sensitivity objectives of upcom- 
ing Ov experiments, a thorough evaluation of the residual 


muon-induced background is imperative. 

To estimate the muon-induced background level, it’s essen- 
tial to understand both the energy and angular distributions 
of underground muons in a given experimental cavern. The 
muon energy was derived using a parametric function [41]: 


dN 
—— = Ce YD (EB + e,(1—e7))-%, 2) 
dE, 


where b = 0.4/km.w.e, Yu = 3.77, and E, = 693 GeV 
were used, as suggested. 
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Fig. 8. (Color online) Energy spectrum of the CJPL muons. Re- 
sults from the parametric function are shown in comparison with 
that from the Monte Carlo simulation. 
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Fig. 9. (Color online)Simulated results of the muon-induced back- 
ground contribution. 


Fig. 8 shows the underground muon spectrum from a 
Geant4-based simulation study and compares it with the Mei- 
Hime parameterization results [42]. The parametric function 
accurately described the Geant4 simulation results. For the 


angular distribution, additional information regarding rock 
overburden is required to determine the flux in a specific di- 
rection. In this study, the zenith and azimuthal angle distri- 
butions of the underground muons were sampled according 
to JNE measurements [43]. The experimental distributions 
can be described using the parametric function discussed in 
Ref. [44]. The total underground muon flux used as the input 
for the simulation was 3.53 x 1071? em~?s~! based on the 
measurement at CJPL Phase II. 

The simulated results of the muon-induced background 
contributions are shown in Fig. 9. We observed that the back- 
ground contribution of muon events, including muon-induced 
secondaries, corresponds to a low level of 2.2 x 1075 ct- 
s/keV/kg/yr in the !°°Mo 0v 66 ROL It was also observed that 
the multiplicity anti-coincidence cut significantly suppressed 
the muon-induced background by an order of magnitude. 


IV. CONCLUSION 


We conducted a simulation study evaluating the back- 
ground shielding for the Jinping !°°Mo-based bolometric 
demonstration experiment aimed at searching for neutrino- 
less double-beta decay. Having established the shielding 
arrangement through an attenuation study, we assessed the 
environmental gamma, neutron, and muon-induced back- 
grounds. Our findings indicate that within the '°°Mo 0vGB 
ROI, gamma rays and neutrons contribute to background lev- 
els of approximately 1.3 x 1073 cts/keV/kg/yr and 1.6 x 
10~° cts/keV/kg/yr, respectively. The underground cosmic 
ray muon’s contribution was minor, only reaching 2.2 x 1075 
cts/keV/kg/yr in the ROI when the anti-coincidence cut was 
applied. Additionally, given the bolometer’s slow data acqui- 
sition rate, potential pile-up events, notably coincident low- 
energy gamma rays, could also add to the discernible exper- 
imental background. It’s worth noting that our simulations 
were done for a demonstration experiment without supple- 
mentary passive background shielding. In CJPL Phase II, a 
meticulously designed combined shielding apparatus (com- 
prising 20 cm of lead and 100 cm of PE) was implemented 
for radioactive background protection. We anticipate a signif- 
icant reduction in the combined contributions from environ- 
mental gamma and neutron backgrounds to levels below 107° 
cts/keV/kg/y, aligning with our objectives for this demon- 
stration experiment. For the 016 experiment, background 
plays a pivotal role, directly influencing experimental sensi- 
tivity. Therefore, curbing background contribution emerges 
as a critical facet of bolometric experimental studies and ad- 
vancements. Regarding the environmental background, pre- 
cise in-situ measurements of both neutron and gamma spec- 
tra, extending into high-energy regions, are paramount for de- 
signing shielding systems tailored to future high-sensitivity 
Ov G6 research endeavors at CJPL 
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